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Abstract 



We have investigated the accuracies of calculations made by omitting the higher 
partial waves of nuclear wave functions and the elementary relativistic terms in 
the hypertriton electroproduction. We found that an accurate calculation would 
still be obtained if we used at least three lowest partial waves with isospin T = 0. 
Furthermore, we found that the omission of the relativistic terms in the elemen- 
^ I tary process amplitude could lead to a large deviation from the full calculation. 

C ^ ■ We also present the cpu-times required to calculate the cross sections. For future 

consideration the use of these lowest partial waves is suggested, since the cal- 
culated cross section deviates only about 0.17 nb/sr 4%), at most, from the 
^ I full calculation, whereas the cpu-time is reduced by a factor of 60. Comparison 

■ of our result with the available experimental data supports these findings. 

' Key words: Meson electroproduction, partial waves, hypernuclei 

' PACS: 13.60.Le, 25.30.Rw, 21.80.+a 



1. Introduction 



An accurate calculation using a simple formalism is naturally desired in all 
KJi \ phenomenological studies of nuclear and particle physics. However, in most 

}_j ■ cases this is difficult to achieve, because an accurate calculation usually does 

I not allow for any extreme approximation. Therefore, an optimal approximation 

which only includes the most important parts of the formalism, without sacri- 
ficing the accuracy of the numerical result, should be obtained. A good example 
is found in the analysis of the hypertriton photo- and electroproduction, i.e., 

7 + ^Hc ^ if + + lU and e + '"^He ^ e' + K+ + , 

given in Refs. [H, 0|- In this analysis the cross section is calculated by means of 
the elementary operator Kaon-Maid W\ sandwiched between two nuclear wave 
functions ('^He and ^H) obtained from the solutions of Faddeev equations using 
modern nuclcon-nucleon and hypcron-nuclcon potentials 0, 0] . The numbers of 
partial waves for the '^He and wave functions are 34 and 16, respectively. 
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Along with the corresponding probabiHties these partial waves are shown in 
Table [TJ In both wave functions the total numbers of supporting points for the 
two-body (p) and the spectator (q) momenta are 34 and 20, respectively. As 
we shall see in the next Section, this could lead to a problem of integration with 
almost two billions grid points. Therefore, it is obviously very important to 
limit the number of participating partial waves or to truncate the elementary 
amplitude in order to simplify the formalism as well as to avoid the unnecessarily 
long cpu-time required to calculate the reaction cross section. At W = 4.04 
GeV it has been shown in Ref. that for the hypertriton photoproduction the 
use of four lowest partial waves {a < 4, see Table [T] for the explanation of a) 
would nicely approximate the full calculation, whereas the use of a < 5 would 
lead to a perfect result. It has been also pointed out that careful inspections 
in a wide range of kinematics should be performed, before we can apply this 
approximation in the hypertriton photo- and electroproduction [H . 

Furthermore, it has been also known that the hypernucleus production cross 
section is sensitive to the elementary amplitude, especially at the forward direc- 
tions, where the two recent experimental data sets from SAPHIR 6] and CLAS 
collaborations show a lack of mutual consistency [S]. As a consequence, hy- 
pertriton production at this kinematics could also shed light on the solution of 
this discrepancy problem. However, the extraction of the information on the 
elementary amplitude from the nuclear cross sections requires a massive fitting 
process, which would become impossible if the cpu-times required to calculate 
these cross sections were extremely long. 

The present analysis is greatly motivated by these facts. Here we shall quan- 
titatively investigate the effects of omitting the higher partial waves of nuclear 
wave functions on the accuracy of the calculation. We shall also compare this 
result with the result of the approximation made by excluding the relativistic 
terms in the elementary amplitude as suggested in Ref. [9[. To this end we 
take the electroproduction process, since photoproduction is only a special case 
of electroproduction. Our motivation is obvious, namely to find the shortest 
cpu-time for which the deviation of the calculated cross section from the full 
calculation is still controllable. 

This paper is organized as follows: In Section [3 wc shall briefly review our 
formalism and state our problem. We shall present and discuss the results of 
our calculations in Section [3l Section |4] will focus on the comparison of our 
results with experimental data. In Section [5] we shall summarize our findings. 



2. Formalism and Statement of the Problem 

The formalism of the hypertriton electroproduction off ^He in an impulse 
approximation has been presented in Ref. [l[ . Except for the difference between 
the initial and final nuclear masses, as well as between the initial proton and 
the final hyperon masses in the elementary operator, the formulas are similar to 
those used in pion electroproduction off ^He [13]. To facilitate the discussion, 
here we will only present the most important part of them. We start with the 
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Table 1: Quantum numbers and probabilities (in %) of the "^He and the hypertriton wave 



functions 4, 5]. 
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corresponding nuclear transition matrix element, which can be written as []| 
(iH| J^'He) = ^/3 J d'p d'q ^^{p,q') J^^ {k,k,,k[) ^s^,{p,q) , (1) 

3 



. K (qic) 




Figure 1: Electroproduction of the hypertriton on a ^He target in an impulse approximation, 
where the virtual photon interacts with only one nucleon inside the '^He. 



where the factor of \/3 on the right hand side of Eq. ([T]) comes from the anti- 
symmetry of the initial state, J** represents the elementary operator, while the 
integrations are taken over the three-body momentum coordinates (see Fig. [T] 
for the explanation of the momenta) 

P = ^ (^2 - fcs) , q = fei , (2) 
and the hyperon momentum in the hypertriton is given by 

q' = ki + ^{k- q^) . (3) 

By expanding the three-body wave functions in Eq. ([1]) in terms of their 
orbital momenta, spins, and isospins, 

*(P:9) = '^(j)a{p,q){LmLSms\Jmj){lmi\ms\jmj){Jmjjmj\\M{) 

Q:,m 

xi"™Jp)i"™,(9)x^.xi \ {T\)\Mt) , (4) 
we can recast the transition matrix element in Eq. ([1]) into the form 

(iHlJ'^I^He) = %/3 ^^(Lmi5ms|Jmj) 

a, a' mm' 

X {LmLSms\J'mji) (^lmi^ms\jmj'j 
X {l'mi,\ms>\j'mj,) [jmjjmj\\M-^ 

^SS' ^msm^, Sto 
X J p^dp d\ (p, q') 0„ (p, q) yS,^, (q') F,^, (q) 

X (i,mv| Jni"^.), (5) 
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where a and a' indicate the partial waves of the initial and final nuclear wave 
functions, respectively, while the notations m = {niL'msrnims'mjmj) and m' = 
{mL'ms'mi>ms'mj'mj') have been introduced for the sake of brevity. 

From Table [T] we can estimate that a full calculation of the four-dimensional 
integrals in Eq. ([5]) using all partial waves could involve integrations with 34 x 
16 X 34 X 20 X 30 X 10 = 110, 976, 000 grid points, where the two last numbers in 
the multiplications (30 x 10) come from the minimum Gauss supporting points 
for the numerically-stable angular integrations [ij. 

It is also important to note that these integrations are performed over all 
components of the transition matrix element in the form of 4 x 4 matrix [ ]ml 
(or equivalently [j^i/]), where 



1 

E 

n— m. 



(1, (Ta;, (Jy, (Tz) 



/ joO jxQ jyO jzQ \ 

JQx Jxx Jyx Jzx 

JOy Jxy Jyy J zy 

\ joz jxz jyz jzz j 



(6) 



As a consequence, the problem of calculating the cross section becomes numer- 
ically more challenging, since it is equivalent to the problem of integration with 
1,775,616,000 grid points. Furthermore, the result of this integration must be 
summed over angular-momentum and spin projections mj ,mji ^ms, and mg 
[indicated by m and m' in Eq. (O]. Fortunately, the selection rule represented 
by the three Kronecker delta functions in Eq. ([5]) along with current conserva- 
tion reduce this number to about 156 millions grid points. Nevertheless, this 
still indicates a time-consuming numerical computation. 

On the other hand, the elementary operator of the elementary process 7i,(fc)-|- 
P{Pp) ^ K^{qK) + A(pa) can be written as 



Fi cr ■ e + F2 cr ■ k €o + F3 cr ■ k k ■ € 

+F4 a ■ k Pp ■ e + F5 a- ■ k pj^ ■ e + — ji^e o- ■ Pp^o + F7 cr ■ p^k ■ e 

+Fs (T ■ PpPp ■ e + Fq(t ■ PpPf^ - e + Fi4(T ■ ea- ■ k(T ■ Pp^ 

+ — \FiQeT ■ p^eo + Fii(T ■ p^k ■ e + Fi2(T ■ p^p ■ e 

+Fi3cr pf^p^ - e + Fi5cr p^a- ■ €(T ■ k^ 

H —(Fiqct ■ p^cr ■ ecr ■ Pp + Fu er ■ p^cr ■ k a ■ p^eo 

£p Ea 

+Fiscr ■ Pf^cr ■ kcr ■ Ppk ■ e + Fig cr ■ p^cr ■ k cr ■ PpPp ■ e 
+ F2Q a- ■ p^a- ■ k cr ■ PpPf^ ■ Xp , 



(7) 
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where the individual amplitudes Fi are given in Refs. [J is the virtual 

photon polarization vector, Sp — Ep + nip, and ea = Ea + tua- 

From Eq. ([7]) it is obvious that the amplitudes Fig — F20 and Fg — F15 
originate from the "small-small" (SS) and "small-big" (SB) terms of the Dirac 
spinors, respectively, whereas the rest correspond to the "big-big" (BB) terms. 
It can be easily shown that the ratio between the "small" and "big" terms is of 



the order v/c (see, e.g., Ref. [12|). In pion photoproduction near threshold it has 
been widely known that only the leading KroU-Ruderman term Fi significantly 
contributes to the process. However, in kaon photo- and electroproduction the 
high threshold energy of the process could certainly change this picture. In 
view of this, it is certainly tempting to neglect the relativistic SS-terms in our 
formalism as well as to investigate the effects of the SB- and other terms on the 
calculated cross sections of the hypertriton electroproduction. 

For the purpose of numerical computation we note that the cross section of 
the hypertriton electroproduction in the cm. system can be written as 

dn^ = dn^ + ''^dn^ + '^ + ^2.l(i + .) ^ cos^^ , (8) 

with e represents the virtual photon polarization, = — (fc^/A;^)e and the 
individual cross sections can be written as 

Jr^-^^f ' (-T,L,TT,andLT), (9) 

with Q!e — e^/47r and Kl = {W'^ — M^^) /2Mb_c- The nuclear structure functions 
Wi are given by 

Wt ^ ^ (VF^^ + wyy) , (10) 

47r 

Wi^ = ^ W°° , (11) 
47r 

Wtt = {W'' - v/yy) , (12) 

47r 

Wist = {W^"" + , (13) 

47r 

where the spin averaged Lorentz tensor W^'^ is related to the nuclear transition 
matrix element given in Eq. ([1]) by 

= ^E(iH|Jn^He)(iH|jn^He)* . (14) 



The kinematic chosen in this analysis is close to that of experimental data , 
because we want to explore the kinematics region that is experimentally acces- 
sible. The angular distribution will be limited in the range of 0° < 6*^ < 31°, in 
which the magnitude of the cross section is sizeable. The total cm. energy is 
limited to 3.5 GeV <W<b GeV, since below this limit the cross section would 
be too small, whereas above this limit the elementary operator would lose its 
predictive power. 
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Figure 2; Effects of tfie higher partial waves on the differential cross section of the hypertriton 
electroproduction off ^He at = —0.35 GeV^ and e = 0.762. The upper-left panel shows the 
result obtained from the full calculation by using all partial waves. Other panels display the 
differences between the full calculation and the calculation by using a < 5 (A5), a < 4 (A4), 
and only s-waves (As). Note that for the sake of visibility a different scale has been used for 
the vertical axis of the upper-left panel. 

3. Results and Discussion 

Figure demonstrates the cross section obtained from the full calculation 
using all partial waves and the deviations from this result if we use a < 5, a < 4, 
and s-waves (a = 2,4). Our experience shows that computing the numerical 
data required by the plot of the differential cross section shown in the upper- 
left panel of Fig. [5] (consisting of 31x31 = 961 points) on a PC with a single 
processor Pentium-4 takes about 11 days (15,344 min). 

From the lower-left panel of Fig. [2] it is obvious that limiting the partial 
waves up to a = 5 yields an accurate approximation, since in general it just 
slightly underestimates the full calculation. The largest discrepancies are found 
at the two cross section peaks at W « 3.75 GeV and 4.10 GeV close to the 
forward angle, i.e., about 0.15 nb/sr (less than 3%). The average deviation over 
these 961 points is only 0.019 nb/sr. A similar behavior is also found if we use 
a < 4, except in this case the calculated cross section slightly overestimates the 
cross section of the full calculation. Here, the average and largest deviations 
are found to be 0.038 nb/sr and 0.17 nb/sr (less than 4%), respectively. Finally, 
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Figure 3: Effects of the relativistic terms on the differential cross section of the hypertriton 
electroproduction off ^ He. The upper- left panel shows the result of a full calculation obtained 
by using all terms in Eq. jTj, i.e., ^1-^20- Other panels display the differences between the 
full calculation and the calculation by using Fi-Fi^ (Agg), F1-F5 (Asb), and only Fi (Ap-^). 
Note that all figures are obtained by using a < 5 and the same kinematics as in Fig. [2] For 
the sake of visibility different scales have been used in vertical axes. 



the largest deviation, almost 1 nb/sr at the top of the highest cross section 
peak, is obtained if we use only s-waves. Since the largest differential cross 
section is around 5 nb/sr, it is obvious that the latter provides a relatively poor 
approximation method for the hypertriton electroproduction. 

Figure [3] demonstrates the effects of excluding the relativistic terms on the 
calculated differential cross section. Since it has just been shown that the use 
of a < 5 can approximate the full calculation with an accuracy up to 3%, we 
believe that it is sufhcient to limit the use of partial waves up to a < 5 at this 
stage. 

The lower-left panel of Fig. [3] demonstrates the effect of the SS-terms exclu- 
sion. We note that the maximum deviation is only about 4% at « 3.75 GeV. 
Therefore, the SS-terms can be safely neglected for this process. 

The shape of the differential cross section obtained by excluding the SB- 
terms is found to be completely different. As a consequence, the deviation from 
the full calculation is quite large (see also Table [5]). The extreme nonrelativistic 
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Table 2: The cpu-time (t) required to compute the numerical data of the cross section plot 
shown at the top panel of Fig. [2] along with the average and maximum deviations from the 
full calculation (indicated by Aav. and Amax. , respectively) for different approximations, i.e., 
using all partial waves (Full), a < 5 (1), a < 4 (2), and only s-waves (3). With a < 5 the 
role of different elementary amplitudes Fi given in Eq. (0 is demonstrated, i.e., using only 
Fi - Fi5 (a), Fi - Fs (b), and Fi (c). 





Full 


1 


2 


3 


a 


b 


c 


T (mill) 


15,344 


511 


294 


46 


498 


487 


485 


Aav. (nb/sr) 




0.02 


0.04 


0.14 


0.02 


0.98 


0.81 


Amax. (nb/sr) 




0.15 


0.17 


0.99 


0.10 


8.03 


4.97 



approximation {Fi only) also yields a large deviation from the full calculation. 
The reason that the deviation in the latter is relatively smaller than in the 
former is that this approximation yields very small cross section. Although 
the analysis of Ref. [9| has been performed at the elementary level and using 
a quite different elementary operator, our result corroborates its finding, i.e., 
both the exclusion of the SB-terms and the extreme nonrelativistic approach 
can not be used as a good approximation in the hypertriton electroproduction. 
Furthermore since the exclusion of the SB-terms was not investigated by the 
authors of Ref. Q , our analysis therefore provides an extension of their finding. 

At this stage, it is also important to consider the cpu-times required to 
make the plots just shown. Along with the corresponding deviations from the 
full calculation, the required cpu-times to calculate the 961 points of the cross 
section by using a single processor Pentium-4 PC are listed in Tabled 

It is interesting to see that the cpu-time is significantly reduced by a factor of 
30 if we limit the partial waves up to a = 5, while the accuracy is still maintained 
up to about 0.15 nb/sr. As a consequence, to obtain the plot shown at the top 
panel of Fig. [3] we need less than 9 hours. If we used the partial waves with 
a < 4, the cpu-time is reduced by a factor of about 60, whereas the maximum 
deviation slightly increases to 0.17 nb/sr. The use of only s-waves substantially 
reduces the cpu-time, i.e., by a factor of 300. The average deviations displayed 
in Table [2] show in general the same behavior. 

In contrast to the omission of the higher partial waves, the exclusion of 
the relativistic terms in Eq. ([7]) does not significantly reduce the cpu-times. 
This is clearly demonstrated in the first line of Table [2l where the cpu-time 
for the calculation with a < 5 decreases from 511 min (full terms) to 485 min 
(only Fi). This result proves that computing the rest 19 Fi amplitudes in 
Eq. ([7]) is much simpler, and therefore much faster, than computing the massive 
integrals given in Eq. ([5]) with the complete partial waves. At the same time the 
maximum deviation in the latter becomes almost 5 nb/sr. This is because in 
the electromagnetic production of kaon the high energy of the process enhances 
the role of all but the Kroll-Ruderman term [see Eq. ([7])], and the extreme 
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Figure 4: (Color online) Comparison between experimental data [Til and the calculation 
using all and specific numbers of partial waves (upper panel) and the approximations made 
by including only certain elementary amplitudes in the elementary operator (lower panel). 



nonrelativistic approach results in a very small cross section. Consequently, 
the largest difference with the full calculation is about 5 nb/sr, i.e., the largest 
differential cross section found with the full calculation. Note that the decrease 
of Amax. from 0.15 nb/sr (full terms) to 0.10 nb/sr (using only Fi — F15) seems 
to be fortuitous, because the average deviations for the two cases are the same 
(0.02 nb/sr). 

To complete our analysis we have also investigated the effects of the higher 
partial waves and the relativistic terms at fc^ = — 1.0 GeV^. However, since the 
obtained cross sections are quite small (0.32 nb/sr, at most), it is very hard to 
draw a quantitative conclusion at this kinematics. 

4. Comparison with experimental data 

Although our primary motivation is to quantitatively study the effects of 
excluding higher partial waves and relativistic terms, it is also imperative to 
compare the results with the available experimental data given in Ref. [Tlj . 
The comparison is displayed in Fig. 3) Unfortunately, due to their large error 
bars, the present experimental data still allow for the s-wave approximation and 
the use of only BB terms in the elementary amplitude. In the latter, it is obvious 
that experimental data at forward angles with about 10% error bars would be 
able to justify the important role of the SB terms Fq—Fi^. In the former, we note 
that the largest deviation of using only s- waves does not appear a.t W = 4.10 
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GeV. Instead, at the forward directions the largest deviation (approximately 
0.99 nh/sT, see Fig. [2]) is found with W ~ 4.20 GeV. Therefore, experimental 
data with the same quality, but atW — 4.20 GeV, would be very useful to check 
the validity of the s-waves approximation. For the sake of numerical accuracy 
and efficient cpu-time we would, however, recommend the use of partial waves 
with a < 4 along with a full elementary amplitude. 

The discrepancy between the calculated cross sections and the experimental 
data point aX 9 — 18.9° requires a special explanation. At this kinematics we 
note that the calculated cross sections are much smaller than those at the for- 
ward direction. This behavior seems to be almost independent of the total cm. 
energy. However, we have also found that the longitudinal part of the cross 
section {daj^/dfl) dominates other parts [see Eq. ([8])] in the whole kinematics 
shown in Fig. [2] and falls off quickly as a function of the 9. On the other hand, 
the angular distribution of the transverse cross section {daT/dil) tends to be 
more flat than the longitudinal one. As a consequence, we may conclude that 
such behavior should originate from the elementary amplitude and not from 
the effect of the nuclear wave functions. Thus, if we believed that this experi- 
mental data point were correct, then we had to reconsider the improvement of 
the elementary operator. As stated in the Introduction, this means that the 
extraction of the elementary information from the hypertriton production cross 
section would be mandatory in the future works. Otherwise, new measurements 
at this kinematics are urgently required. 

5. Conclusion 

In conclusion, we have investigated the effects of higher partial waves and 
relativistic terms on the accuracy of the calculated differential cross sections of 
the hypertriton electroproduction. We have shown that an accurate calculation, 
with a maximum deviation of less than 4%, could still be obtained if we used 
the three lowest partial waves with isospin zero (i.e., using a < 4, since the 
selection rule excludes the a = 1 component). In this case, the cpu-time for 
calculating differential cross sections is reduced by a factor of about 60. The 
exclusion of certain elementary amplitudes Fi has a tiny impact on the cpu- 
time, but a big impact on the accuracy of the calculation. In view of this, 
for future consideration we suggest the use of partial waves with a < 4 with 
a full elementary amplitudes Fi. Comparison of the results with the available 
experimental data supports our finding. New measurement of the hypertriton 
electroproduction with about 10% error bars would be very useful to clarify the 
validity of these approximations as well as the importance of the relativistic 
terms. 
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